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ABSTRACT 



The spin evolution of the compact neutron core in a Thorne-Zytkow Object 
(TZO) is investigated to explore the origin of extremely long period X-ray source. 
It is found that the outflow would effectively take away angular momentum from 
the core when radiation pressure dominates the accretion process. Thus the 
compact core could quickly spin-down to the co-rotation period (e.g. several 
hours) within the massive envelope, in about 10 3 — 10 4 years. The compact core 
could become an extremely long period compact star if the envelope is disrupted 
by some powerful bursts or exhausted via the stellar wind. The 6.67-hour periodic 
modulation of the central compact object (CCO) in supernova remnant RCW 103 
could be naturally understood as the descendant of a TZO. 

Subject headings: ISM: supernova remnants - stars: evolution - stars: individual 
(1E161348-5055) - stars: neutron 
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Introduction 



The central compact object (CCO) in supernova remnant RCW 103, 1E16 1348-5055 
(hereafter 1E1613), is a very strange X-ray point source. It has a X-ray luminosity in the 
range of ~ 10 33 — 10 35 erg/s, and has no identified radio, infrared or optical counterpart. 
In 2005, a periodic mod ulation of 6.67 hours was detected by the EPIC MOS cameras 



onboard XMM-Newton ( De Luca et al. 



20061 ). and any other periodicities with P > 12 ms 



are excluded with hig h confidence. However, RCW 103 is a supernova remnant of only 



1984 



Carter et al 



19971 )). which means that the CCO 



about 2000 years old (INugent et al.l 
E1613 is very young. These features make 1E1613 a unique source among all the CCOs 



(IDe Luca 



2008|). 



The 6.67-hour periodicity is normal if it is the orbital period of a binary. But it is 
difficult to explain how the companion survive from the supernova explosion and escape 
from the multi-wavelength observations, and why the spin period of 1E1613 has not been 
detected if 6.67 hour is the orbital period. If 6.67 hour is the spin period of 1E1613, it 
would be the longest one of all is olated pulsars. How c ould an isolated pulsar spin-down to 



such a long period in 2000 years? IDe Luca et al. 



(120061 ) proposed a scheme that a magnetar 
with an initial period of 0.3 s and a surface magnetic field of 5 x 10 15 G surrounded by a 
disk of 3 x 1O _5 M may spin down to a period of 6.67 hour in 2000 years, owning to the 
extrem ely large propeller torque. T his initial period is too large according to the magnetar 



theory (IDuncan fc Thompson 



19921 ). and its current period is very different from the range 
of 2-12 s for the other magnetars detected as ano malous X-ray pulsa rs (AXPs) and soft 



gamma-ray repeaters (SGRs). The AXP 0142+61 (IWang et al.l 120061 ) has a period of only 
8.7 s. We noticed that even in binaries, extremely long periods have been detected from 
X-ray sources in binaries, such as the ^2.7 hour period X-r ay source in 2S 0114+65 



(IFinley. Belloni. fc Cassinelli 



1992 



Li fc van den Heuvel 



19991 ) and the ~ 1.5 hour period 



-4- 



X-ray source in 4U 2206+54 flWaneJl2009h . 



On the other hand, the predicted star s wit h compa c t neu tron cores remain mystery 



which was first proposed by iGamowl (jl937l ) and lLandaul (jl938f l to solve the stellar energy 
problem. The idea was aband oned when t h e stellar energy is r ealized to b e generated 



from thermonuclear reactions (IBethe 



19391). 



Thorne fc Zytkowl (11975 



19771 ) discussed the 



possible equilibrium states of a massive star with a degenerate neutron core, which was 
called Thorne-Zytkow Object (TZO). A TZO could be produced in at least three ways: 1) 
The expansion of the normal star in a binary drags the companion neutron star (NS) and 



gradually merges (Taam et a 



star (MS star) (IBenz k Hills 



1978); 2) A collision of a NS with a massive main-sequence 



19921 ); and 3) A supernova (SN) explosion gives the young NS 



an appropriate kick ve locity and a direction whic h make it embed in its binary companion 



( [Leonard et al. 



19941 ). 



Podsiadlowski et al 



(119951) estimated a total number of 20-200 TZOs 
existing in the Galaxy given the characteristic lifetime of TZO as being 10 5 — 10 6 yr. 



The TZOs usuall y manifest as red supergiants or red giants (IThorne fc Zytkow 



1977 



Eicheta] 



19891 ). But they can be distinguished from the ordinary ones since 



they have extraordinary high abundances of lithium and rapid proton process elements 



( IThorne fc Zytkowl 1 19771 ). The final fates of TZOs are not entirely clear. It is generally 



believed that a neutron core accreted enough material and finally transform to a black 
hole (BH), and that the massive envelope may be dissipated and form a disk. Here we 
discuss another possibility. If a powerful burst occurs before the neutron core transform to 
a BH, the burst can disrupts the envelope, the neutron core may looks like a CCO, and the 
envelope may mix with the original SNR which produced the neutron core or form a new 
SNR. Usually, the spin energy of the neutron core is much smaller than its kinetic energy 
and the accretion energy, thus has little influence to the evolution of a TZO. However, if the 
core could be observable after the envelope is disrupted, the spin behavior of the core could 
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be very important. We proposed that 1E1613 and RCW 103 are the TZO descendant. We 
will deliberate how the compact core in a TZO spin down to a long period, and analyze 
different fates of the compact cores in TZOs. There may be two possible approaches for a 
TZO transforming to 1E1613 and RCW 103. 

Case A. About 2000 years ago, the more massive star of a binary experienced a 
supernova explosion, which produced RCW 103 and a compact star. By some possibility 
the compact star got an appropriate kic k velocity and direct ion to be embedded in the 



companion MS star and formed a TZO flLeonard et a 



as a red giant or red supergiant afterwards (jEich et al 



1994). The TZO would manifest 



19891 ). which has angular frequency 



much smaller than its compact core. Thus, the core is then spun down to the co-rotational 
period (6.67 hour) by the envelop efficiently, as discussed below. And then a powerful 
burst occurred and destroyed the envelope (e.g., due to a phase transition of its crust 



(jCheng et al 



19981 ) . or other processes such as thermal nuclear explosion). As a result, the 
envelope mixed together with the original SNR RCW 103, and the core became observable 
as a central compact object, i.e. 1E1613. The accretion from the residual envelope or a 
fall-back disk reproduces the observational X-ray emission. 



Case B. More than 2000 years ago, a neutron star in a 



its companion star or was swallowed by it (iBenz fc Hills 



Dinary was embedded in 



19921 ) and formed a TZO. The 



core spin down quickly to the co-rotational period, 6.67 hours. About 2000 years ago, a 



global thermal nuclear exp 



osion took p 


ace o 


(Ouved et al. 


2002) 



ace o r the whole neutron star occurred a phase 



20021 ). The powerful energy emission destroyed the 



TZO envelope and formed an analogous supernova remnant, i.e. RCW 103. The central 
compact object 1E1613 accretes from the residual envelope or a fall-back disk and emits 
the observational X-ray. 



In both cases above, the central compact object would has very small proper motion 
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velocity, which is consistent with the observation of 1E1613. 



2. Spin evolution of the neutron core in a TZO 



Podsiadlowski et al 



(119951 ) estimated the spin evolution of a neutron core in a TZO 
via the total angular momentum carried by the accretion material. They found a slowly 
spin-down process during the steady-burning phase where the accretion has an Eddington 
rate, and a quickly spin-up process in the neutrino runaway case where the accretion 
is not Eddington-limited which could spin-up the core to a typical period of ~ 10 ms. 
Nevertheless, they did not consider the interaction between the neutron core magnetosphere 
and the TZO envelope, and also ignored the outflow of the convectio n envelope. The 



conve ction passes through the whole envelope out to the photosphere ( IThorne fe Zytkow 



19771 ). The outflow carries away the heat energy as well as the angular momentum, which 
would make the neutron core spin-down effectively. We now deliberate the spin evolution 
of the neutron core by considering both the interaction and the outflow. 



2.1. The model 

During the steady-burning phase of a TZO, some of the envelope material is accreted 
to the neutron core surface. The radiation pressure opposes the gravitational force and 
forms an equilibrium. The difference between the inflow rate (M in ) and outflow rate (M out ) 
determines the rest accretion rate (M acc ), i.e. 

Mace = M in - M out . (1) 

In a spherical symmetry accretion case, M acc approximates the Eddington rate, i.e. 

M acc « 10 18 i? 6 g/s, (2) 
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where Rq = -R*/10 6 cm, where R* is the compact star radius. 

To estimate the inflow rate we consider the accretion in the steady medium, whose 
accretion rate is ~ nRlc s p, where R a is the accretion radius, c s and p are the sound velocity 
and density at R a . The sum of the gravitational potential energy and the kinetic energy of 
a particle at R a is zero. When r < R a , the kinetic energy of the particle cannot overcome 



the gravitation, thus the particle is cap tured by the core. 



by a envelope with a structure gave in 



or a neutron core surrounded 



Thorne Zytkowl (119771 ). the accretion rate in a 



steady medium is about 10 21 g/s, much larger than the rest accretion rate. The inflow rate 
should be equal the accretion rate in a steady medium at R a , if the radiation pressure is 
insignificant. When considering radiation pressure, we take 



where < rj < 1. If we set 



we get 



and 



M in (i?a) = V^R^CsP, (3) 



M in (i?*) = M acc , (4) 



Mout(#*) = 



M ont (R a ) = M in (R a ) - M acc . 

In the range of i£* < r < R a , the outflow can be different for different TZOs, we take an 
exponential form of: 

^out(r) = ( 5 - R * ) a (M in (R a ) - M acc ), (5) 
where the parameter a ~ 1. 

The above equations provide the distribution of the inflow and outflow inside R a . 
Their radial velocities are determined by the gravitational force, radiation pressure and 
the distribution. Their rotational velocities, which are related to the angular momentum, 
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should be determined by the interaction between them and the magnetosphere of the 
neutron core. 

The accretion in a TZO is different from the traditi onal accretion of a NS. The envelope 



has extremely slow rotation (IPodsiadlowski et al. 



19951 ) thus a disk could not form around 



the neutron core. The interact ion between the m agnetosphere and accretion flow would not 



19731 ) since the structure of the accretion region 



confined by the Alfven radius (ILamb et al 
is radiation-pressure dominant. 

The magnetosphere of the neutron core is insignificant in the radial direction since the 
convection flows are radiation-pressure dominated. Nevertheless, it would dominate the 
rotational movement of the flows because the gravitational force and radiation pressure do 
not work in this direction. If the flows rotate so fast, that the centrifugal force can compete 
with the gravitational force, the structure of accretion region would be affected. If the 
centrifugal force is much smaller than the gravitational force, the rotational movement and 
radial movement should be relatively independent. Therefore, we can seperatedly study the 
rotational interaction between the magnetosphere and the flows while their distributions 
are determined by the radial direction forces. 

The rotational velocities of the flows are shown in Figure 1. The flows have a potential 
to co-rotate with the magnetosphere with a co-rotational velocity 

v co = rQ (6) 

because of the magnetic freezing effect, where Q is the angular velocity of the compact core. 
But the magnetic field can draw the flows to a velocity at most with kinetic energy density 
equal to the magnetic energy density, i.e. 

v dr = (3B(^pY l/ \ (7) 
here < /3 < 1, and B is the magnetic field strength. Thus the real rotational velocity 
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should be the smaller one of v co and v^r, i.e. 

Vtf, = min(v co ,v dr ). (8) 

As shown in figure 1, the flows would be drawn by the magnetosphere with Vdr when they 
are far from the core, and they would co-rotate with the magnetosphere when they are near 
to the core. At a radius R e , 

v co (R e ) = V d r(R e ), 

we obtain 

R e = (5B(ATxp)- 1/2 n-\ (9) 
Note that B and p are the functions of radius r. 

When a particle accesses the accretion radius R a , it will be drawn by the magnetosphere 
and rotate around the core with velocity v dr . If it passes through the equal velocity radius 
R e , it would co-rotate with the core. Therefore, the particle should be accelerated in the 
magnetosphere and get angular momentum from the core. If a particle finally hits to the 
core surface it would return the angular momentum to the core, as well as the initial 
angular momentum when it was in the envelope. But if it goes out it takes away angular 
momentum from the core, and helps the core to spin down. 

If a particle goes out from R e < r < R a , the angular momentum it takes away from the 
core is 

j = mrv dr , (10) 

where m is the particle mass. If a particle goes out from R* < r < R e , it would take away 
an angular momentum of 

j = mR c v co (R e ). (11) 

Note this value is calculated at R c not at r, because the particle would co-rotate with the 
core when going out from r to R e . When the particles are going out from R c < r < i? a to 
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R a , they could also get a bit more angular momentum from the core because the field lines 
are sweeping them. We ignore this additional amount of the angular momentum, because 
it is the same order as the initial value and can be degenerated with the model parameters 
(e.g. B and (3). 

If we assume the envelope spins with a angular velocity of fi a at R a , we get the angular 
momentum change rate of the neutron core as 

m = j = i in - j out , (12) 

where 

i in = M iQ (R,)^Rl (13) 

is the angular momentum change rate from inflow and J out the angular momentum change 
rate from outflow. From equations (10) and (11) we know that the angular momentums 
took away by the particles go out from the range of r < R c and r > R e are different, thus 
J out should be the sum of the contributions in the two ranges, i.e. 

Jout = Jr<_R c + Jr>R c , (14) 

where 

j r<Rc = WM out v co (R e )R e dr = r° a ^ M ff (r - R^Qdr, (15) 

and 

j r > Rc = VM out v dr rdr = T a ?zT%t ( r ~ R^m^pV^dr. (16) 

In the above equations, B is defined by the structure of the magnetosphere, p defined 
by the structure of the envelope. Therefore, the spin evolution of the compact core in a 
TZO could be numerically obtained from equations (2), (3), (5), (9), and (12)-(16), with a 
set of given parameters. 



- 11 - 



2.2. The results 

In the flowing calculations, we took a dipole magnetic field for the compact core, 
where the magnet i c mom ent fi = Br 3 . We adopted the envelope structure of TZO given by 



Thorne fc Zvtkowl (Il977l ). with P a « r D _ m = 10 9 cm, = 10 6 cm, and p = 10 12 r~ 2 g/cm 3 
when supposed the exponential density distribution. We took the envelope at R a spin 
period P a = 27r/fi a = 6.67 hours (the modulation period of 1E1613) although the results 
are not sensitive to P a . 

If the core spins too fast (e.g. fli in Fig. 1) the centrifugal force may exceed the 
gravitational force at a radius where the model should be invalid. When the core spins so 
fast that the light speed radius (t"l = c/fl, where c is the speed of light) is smaller than the 
accretion radius (i.e. < R a ), the upper limit of integral in equation (16) should change 
from P a to t"l, since the interaction between particles and magnetosphere is negligible when 
r > rjj. On the other hand, when the core spins too slowly (e.g. Q 3 in Fig. 1) where 
R e > P a , equation (16) does not work and the upper limit of integral in equation (15) 
should change from R e to P a . 

Figure 2 shows the spin evolution of a compact core in a TZO with different core 
parameters and same model parameters, where Pq is the initial period of the core, and 
^30 = /i/10 30 G cm 3 . It shows that a compact core can spin down to the envelope spin 
period P a in about 10 3 — 10 4 years, then co-rotate with the envelope. A core with stronger 
magnetic field could spin down faster than a weaker one. The spin-down timescale is 
approximately in proportion to the initial period. 

Figure 3 shows the spin evolution of a compact core in a TZO with same core 
parameters and different model parameters, where Pq = 0.2 and /130 = 1. It shows that a 
core with smaller a, larger and r] spin down more effectively than the one with a larger 
a, smaller j3 and 77. The spin-down timescale is approximately in inverse proportion to the 
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accretion fraction 77. 

When the core has spun down it co-rotate with the TZO envelope during the 
steady-burning phase. If the neutrino runaway phase occurs, it could be spun up by 
the huge rest accretion flow. If not the neutrino runaway phase does not occur the core 
co-rotate until something else happens. 



3. The fate of the compact core in a TZO 

The compact cores in the TZOs with different initial conditions could have different 
fates. The relation between the initial condition and final fate is not clear. The core could 
transform to a BH if it accretes enough material or becomes a long period isolated compact 
star when the envelope vanished. 

If the core evolved to a BH, its spin information should transform to the specific 



angular momentum of the BH, which may 



accreting from an disk around ( iZhang et al. 



je te sted from the X-ray properties when 



19971 ) . Usually, the BH would spin fast since it 



has accreted enough material with angular momentum. 

It could evolve into an isolated compact star when the envelope vanishes in a relaxed 
or an abrupt way. A TZO usually has strong stellar wind, thereby it may lose a large part 
of the envelope material. The core might be naked when most of the envelope has been 
lost, and the core could finally evolve into an isolated neutron star with an extremely long 
period, which might be also detected as an X-ray source due to remnants thermal energy or 
accretion from the residual envelope. 

Three possible processes might destroy an envelope in an abrupt way. 1) When 
non-steady burning occurs, it would disturb the stable evolution. As its energy is not large 
enough it may not disrupt the envelope, but a series of nuclear burst could speed up the 
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progress. 2) If the strange quark matter is the real base-state of hadron (jWitterull984j ). 
the phase transition might happen in extreme conditions. The hadron material transfers 
to quark material, and the process releases a l arge amount of energy which is enough to 



disrupt the TZO envelope. 



Ouyed et al. 



( 120021 ) proposed that a neutron star could transfer 



to a quark star, which releases an extremely large amount of energy, and appears as a 
quark-nova. If this happens after the neutron core has spun down, the envelope can be 
disrupted. The supernova explosion leaves the quark star as a central compact star with 
long period. 3) When the core ha s accreted enough m aterial, the hadron crust might break 



19981 ). If the transforming crust is larger 



and transform to quark material (jCheng et al. 
than ~ 1O" 5 M , the energy could be enough to disrupt the envelope, which may also 
result in a long period isolated pulsar surrounded by diffuse medium rather than supernova 
remnant. If the former supernova remnant has not disappeared, the disrupted envelope 
would mix with it. 

The compact degenerated star, including BH, should has a small proper motion velocity 
since it was generated inside a TZO. 



4. Discussion and conclusions 

The mysterious TZOs should exist in our Galaxy since a neutron star has various 
possibility to enter into a MS star. They could be verified by direct observations of a TZO 
or indirect observation of a descendant. Unfortunately, no direct observation is available. To 
indirectly verify the existence of the TZOs we need to understand what their descendants 
could be. In this paper, we studied the spin evolution of the compact core inside a TZO, 
and concluded that the core can quickly spin down to the co-rotation period in the inner 
part of the envelope. A TZO could evolve into a long period isolated X-ray compact star, 
surrounded by a supernova or quark-nova remnant. The observations of 1E1613 in RCW 
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103 are in accord with this scheme in the long periodicity, small proper motion, and other 
observed X-ray emission features. 

We have calculated the spin evolution of the compact core in the possible equilibrium 
states of a massive star with a degenerate neutron core (i.e. the TZOs). The compact 
core of a TZO would spin-down quickly as braked by the interaction between the envelope 
and magnetosphere. The core could manifest as a long period isolated compact star if an 
energetic burst destroys the envelope or it is faded away with the stellar wind. The CCO 
in RCW 103, i.e. 1E1613, might be a descendant of a disrupted TZO. A few compact stars 
possess an extremely long period. Most of them may not be observable because they do not 
have enough accretion from surroundings to produce strong X-ray emission. 

The compact stars could also be one of the possible origin of the extremely long period 
objects in binary X-ray sources, if they were captured by other stars. In this scenario the 
extremely stron g magnetic field magnetars a nd their magnetic field decay are not necessary, 



Li fc van den 



as discussed by 
source in 2S 0114+65 and by 
2206+54. 



Heuvell (|1999h for the origin of the ~ 2.7 hour period X-ray 



Wangi (120091 ) for the ~ 1.5 hour period X-ray source in 4U 
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Fig. 1. — The rotational velocity of a particle in accretion region. is the compact neutron 
core radius, R e the equal velocity radius, and _R a the accretion radius where accretion begins. 
When a particle accesses the accretion radius R a , it is drawn by the magnetosphere and 
rotates around the core with velocity t>d r . If it could passes through the radius R e , it will 
co-rotate with the core with velocity v co . When a particle goes out from the range of 
R e < t < Ra, or R* < r < R e , it could not be decelerated, thus would take away angular 
momentum with an amount given in equations (10) and (11), respectively. 
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Fig. 2. — The spin evolution of the compact core in a TZO with different core parameters 
and same model parameters, where a — 1, /3 — 1 and rj = 0.5. We take I— 1 30 gcm 2 for the 
compact core and Pe = 6.67 hours for the envelope. 
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Fig. 3. — The spin evolution of the compact core in a TZO with same core parameters and 
different model parameters, where fi^ = 1 and Pq = 0.2 s. We take /= 10 30 gcm 2 for the 
compact core and Pe = 6.67 hours for the envelope. 



